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1. Introduction 

Colloidal dispersions can be driven into stationary non-equilibrium states by shearing. 
Their properties are important for the handling of dense colloidal dispersions, but yet 
not well understood from fundamental starting points [1]. A prominent and universally 
observed effect is shear thinning, that the viscosity of the solution decreases by orders 
of magnitude upon increasing the shear rate. A widely used many-body model of 
colloidal particles under shear is provided by the Smoluchowski equation [2] , a special 
form of a Fokker-Planck equation [3], which, however, has yielded to exact solution 
only at low particle densities [4]. There it exhibits weak shear thinning followed by 
shear thickening. The model supposes the existence of a given solvent velocity flow 
profile that depends linearly on distance along one direction, as has been observed 
in polydisperse dense fluid and glassy colloidal systems [5,6], and thus it docs not 
include changes of the solvent velocity field induced by the particle motion. (Nor does 
it allow for 'shear banding' or other symmetry breaking phenomena.) An additional 
simplification of the model can be performed when solvent induced interactions 
('hydrodynamic interactions') are neglected, so that the model effectively describes 
interacting Brownian particles in a constant shear flow. Gratifyingly, shear thinning 
has been observed in simulations of this system [7] , where also homogeneous states in 
a linear flow profile were recorded. The model thus contains flow curves (viz. curves 
of stress versus shear rate) in qualitative agreement with typical experimental data 
of systems close to glassy arrest [8] at not too high shear rates (where hydrodynamic 
interactions presumably dominate [1]). 

Recently we presented a mode-coupling approach which leads to a consistent 
and (in principle) parameter-free, quantitative, albeit approximate description of 
the stationary sheared state at high particle concentrations or strong interactions 
[9], It explains the behaviour of dense dispersions under shear from considering the 
competition of local caging of particles [10], which causes slow structural relaxations, 
with shear advection of fluctuations [11], which speeds up the decay. As an important 
concept it uses integration of the transient dynamics in order to gain insights into 
the stationary state presumed to be reached at late times. A number of rather 
universal predictions of the approach have already been obtained [12], and are in 
qualitative agreement [13] with, e.g., the mentioned computer simulations [7]. In 
this contribution we present details of the approach, starting from the many-body 
Smoluchowski equation with shear and setting up the frame for integrating through 
the transient. Formally exact expressions for stationary averages, correlation functions 
and susceptibilities, and for transient density correlators are presented. These exact 
results lay the foundations for our approach to sheared colloids, whose approximations, 
outlined in Ref. [9], will be given in detail in a companion paper. In particular the 
results derived below create a framework within which to make mode-coupling-type 
approximations for sheared colloids, without invoking the equilibrium form of the 
fluctuation-dissipation theorem (which cannot reliably be used under shear). 

It is well known in many physical situations that the same standard approximation 
(e.g., factorising an average) gives different results when applied to two formulations 
of a problem that would, if treated exactly, be equivalent. A careful choice is then 
required, and the work reported here can be thought of as 'preparing the best ground' 
for a judicious mode-coupling approach to sheared colloids. 

The derived generalized Green-Kubo relations and generalized Zwanzig-Mori 
equations, which may also be of interest on their own, are useful to describe the non- 
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equilibrium steady state because they enable one to connect the stationary distribution 
function to the transient dynamics evaluated with equilibrium averaging. This strategy 
was followed in mode-coupling calculations of the non-linear viscosity of simple liquids 
under shear, where the slow relaxation of 'long time tails' leads to non-analytic 
dependences on shear-rate [14] . Transient fluctuation functions were also successfully 
used in some of the simulation studies of this problem, and the connection to the 
theoretical approach was shown explicitly [15]. Here shearing cuts off the anomalous 
long-time dynamics present in the quiescent fluids, and subtle (but small) corrections 
to the viscosity arise. 

Recent mode coupling theory approaches to simple liquids close to glassy arrest 
by Miyazaki and Reichman [16,17], and to violations of the fluctuation-dissipation 
theorem in Brownian particle systems by Szamel [18], follow a somewhat different 
approach from ours. There, time-dependent correlation functions for fluctuations 
around the sheared steady state are obtained, as are susceptibilities describing the 
response of the state in Ref. [18]. Thus, in the spirit of the mode coupling theory of 
quiescent systems [10], structural quantities of the stationary state, which now depend 
on shear rate, enter as input into the equations describing the dynamics. Importantly, 
Miyazaki and Reichman find that shear advection of density fluctuations speeds up 
the structural relaxation, which would become excessively slow close to glassy arrest. 

Our approach, as sketched in Ref. [9], uses generalized Green-Kubo relations 
to access the stationary distribution function, in order to allow for its (possibly) 
non-analytic dependence on shear rate. We expect nonanalyticities to arise when we 
consider the rheological properties of the system close to solidification into a colloidal 
glass, because the quiescent dynamics becomes non-ergodic at a glass transition 
described by the idealized mode coupling theory [10]. Therefore in this contribution, 
(i) in Sect. 3 stationary averages are reformulated so that the transient dynamics 
enters; in the companion paper we will present in detail mode coupling approximations 
so that the transient dynamics is described by transient density fluctuation functions. 
Then (ii) in Sect. ^, the equations of motion of transient density functions are 
reformulated in such a way that dense systems can be described, where particle 
interactions lead to large memory-effects. In Sect. 2 the model is defined and some 
properties discussed, while Sect. 5 concludes with a short outlook to the companion 
paper. 

2. Steady state properties 

2.1. Microscopic starting point 

The system considered consists of N spherical particles (diameter d) dispersed in a 
volume V of solvent with imposed flow profile v(r) = k- r, where for simple shear 
with velocity along the x-axis and its gradient along the y-axis, the shear rate tensor 
is K = 7 icy (viz. K a p = jS ax 5/3 y ). The effect of the shear rate 7 on the particle 
dynamics is measured by the Peclct number [1], Pco = jd 2 /Do, formed with the 
(bare) diffusion coefficient D of a single particle. Dimensionless units are obtained 
by setting d = Dq = kgT = 1, whereupon Peo = 7. The evolution of the distribution 
function ^(T) of the particle positions, r^, i = 1, . . . , N (abbreviated into T = {r,}), 
under internal forces F; = — diU(T) (with the total interaction potential U) and 
shearing, but neglecting hydrodynamic interactions, is given by the Smoluchowski 
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equation [1,2]: 

fl t tt(r ) i) = fi(r)*(r J t) ) 

n =n e + sn = J2 0* ■ (o< - f* - «■ r<) ; (i) 

i 

here f2 e = e?i • (9, — Fj) abbreviates the Smoluchowski Operator (SO) without 
shear. In the following, operators act on everything to the right, if not marked 
differently by bracketing. The conditional probability, for the system to evolve from 
state point V at time t' to V at the later time t, denoted by P(Tt\T't'), also is 
determined from Q: 

d t P(Tt\T't') = fi(r) P(Tt\T't') , (2) 

with the initial condition P(Tt\T't) = S(T — V), that both state points coincide at the 
same time. 

There exist two special time-independent distribution functions; the equilibrium 
one, \l/ e , and the stationary one, which satisfy: 

O e * e = , 0* s = . (3) 

The equilibrium one is determined from the total internal interaction energy U via 
the Boltzmann weight, ^(T) oc e~ ut - r \ but the stationary distribution function ^ s is 
unknown. Averages with \I> e will be abbreviated by (...) = J dT^f e (T) . . ., while *S? S 
determines steady state averages, denoted by (. ..) (-v) = JdT^ s (T)..., 

At finite shearing, steady-state averages /, time-dependent correlation Cf g (t) 
and time-independent structure functions Sf g , for fluctuations Sf = f—{f)™ around 
the steady state, and response susceptibilities Xfg(t)> are the central objects of interest: 



f(i) = (f) W = j dT *,(T) f(T) 

C fg (t) = j dT J dY' W 2 {Yt + t l ,Y't') 6f*(T') 6g(T) 

= (Sf* e nU Sg) {A<) 
Sfg = C fg (t = 0) = (Sf* Sg)W 

i 

The calculation of the fluctuation functions involves the joint probability distribution, 
W2(Tt, T't'), that the system is at point (T,t) after it was in a stationary state at 
(r',<'); it is given using the conditional probability that is the solution of Eq. (2): 

W 2 (Tt,T't') = P(Tt\T't') * fl (r') = e n ( r )(*-*') 5(T - V) V S (T') . 

Exchanging the order of times, t' > t, it obeys W 2 (Tt, T't') = W 2 (T't' 1 Tt). The adjoint 
of the SO arose in the fluctuation function C(t) from partial integrations: 

n+ = 5^(a i + p i + r i -K r )-a i) (5) 

i 

where surface contributions are neglected, throughout, for the considered infinite 
system (V — > oo). The susceptibility Xfg(t) describes the linear change of the 
expectation value of variable g: 

A.g(7) (t) = (g)^ = /* dt> Xfg (t - f) h e {t') + 0(hl) , (6) 
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upon application of an external field h e (t) that couples to the variable /* in the 
potential energy; viz. when the potential energy U is perturbed to: 

U^U-f*(T)h e (t). (7) 

The standard calculations [3] leading from Eq. (7) to Eq. (4) are sketched in the 
Appendix. 

Without applied shear the SO 0| is an Hermitian operator with respect to 
equilibrium averaging [2] 

((nir)fl> = (mu> = -E<^-^>' ( 8 ) 

i 

and (as seen from specialising to / = g) possesses a negative semi-definite spectrum. 
But with shear fit cannot be brought into an Hermitian form [3]; see Sect. 2.2.1 
below. The action of fi on the equilibrium distribution function = 5£lty e will 
become important later on and allows one to define the stress tensor: 

5fi \& e = - 9j k- Ti * e = - V" (Fj>K- r,; + Trace{/c}) * e 

i i 

= Tracc{/« cr} * e = 7 a xy * e , (9) 
with a a p the zero-wavevector limit of the potential part of the stress tensor: 

^ = -E(^+^rf)- (10) 

i 

The specific form of k for sheared systems was used in the last equality of Eq. (9) 
only. 



2.2. Basic properties 

Some well known properties of solutions of Fokker-Planck equations [3] shall be 
collected which bear relevance to the discussion of sheared colloidal dispersions. 



2.2.1. Eigenfunctions expansions: The Smoluchowksi equation of Eq. (1) may be 
viewed as a continuity equation in phase space [3], dt^(T, t) + ^ di ■ J,(r, t) = 0, 
where the probability current equals: 

j<(r,t) = (Fi + K-Ti-di)v(r,t) . (11) 

Stationarity implies J^- di ■ Jj(r, t — > 00) — 0. But only if the current vanishes, J; = 0, 
in the steady state, one can show that the SO fi is related to an Hermitian operator. 
Necessary condition for J; = are the 'potential conditions': 

9f (if + (kt/) = df (F? + (K-r.n (12) 

While the potential conditions hold in equilibrium, they are violated under shear, 
because 9(y • K-Ti))/dxi ^ 9(x- K-r,)) / dyi. Then, if an expansion in eigenfunctions of 
the SO exists, it will have the following properties: The conditional probability from 
Eq. (2) takes the form 

p(rt\r't') = v»(r) ^„(r') e - A - <*-*'> . (13) 

n 

where the eigenvalues satisfy 5RA„ > 0, and the sets of eigenfunctions 

= -A n fn and flT (p n = -\* (p n ; (14) 
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arc bi orthogonal, viz.: 

dT(p n (T) <p m (r) = 5 nm . (15) 



Yet, because no further connection between the sets of eigenfunctions exists in general, 
important properties of equilibrium fluctuations cannot be expected under shear: for 
example the autocorrelation functions Cff(t) of Eq. (4) can fail to be of positive type 
and may exhibit negative frequency spectra [19]. 

2.2.2. Fluctuation Dissipation Theorem: In the case without shear, where f2 e is 
the SO and averages are performed with the equilibrium distribution function a 
simple relation exists between the fluctuation function C^(t) and the susceptibility 
Xfg(t)- By partial integration, and recalling that 9j\& e = F^ e , one finds: 



*/2(*)= - W nU alt 8g) 



= -d t C%\t). (16) 

The expected fluctuation dissipation theorem (FDT) connects response and fluctuation 
function. On the other hand, with shear the susceptibility Xfg(t) is connected to 
a fluctuation function of a variable /, which can only be found if the stationary 
distribution function is known [3, 18]; it satisfies: 

v &i ■ = o ;* * s . (it) 

i 

and the FDT then states [3] : 

Xf B (t) = -(Sf* ^ e nU 6g)M = -dt C fg (t) . (18) 
This appears not to be as useful as Eq. (16). 

2.2.3. Aspects of translational invariance: Homogeneous amorphous systems 
shall be studied so that by assumption the equilibrium distribution function ty e 
is translationally invariant and isotropic. As shown in Sect. 3.1, the steady- 
state distribution function with shear, then also is translationally invariant, 
assuming that no spontaneous symmetry breaking takes place, but anisotropic. 
Appreciable simplifications follow for steady-state quantities of wavevcctor-dcpcndcnt 
fluctuations: 

/ q (r,i)^ ft £ X{(T)e^, (19) 

i 

where e.g. Xf = 1 describes density fluctuations £? q (i), while Xf = S a p + 
(l/2)Y^(rf — r")du(\ri — Tj\)/drf gives the stress tensor element ^^(q) for 
interactions described by the pair-potential u. Translational invariance in an infinite 
sheared system dictates that averages are independent of identical shifts of all particle 
positions, r — > V where = r, + a for all i. Under such a shift the SO becomes 

rjt(r) = fit(r") -a-« T - p , with p = ^9,. (20) 

i 

Thus a fluctuation of a variable which depends on particle separations only, viz. 
X{ (r) = X{ (V) so that PX/(r) = holds, transforms to 

/ q (r,t) = e -'^- K *)- a / q (r',i). (21) 
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As the integral over phase space must agree for either integration variables T or 
r', steady-state averages from Eq. (4) can be non-vanishing for zero wavevector only: 

/o(7) <S q ,o = \ (/ q (i)) W • (22) 

In the following the index will often be suppressed in e.g. the average density 
g = N/V and the shear stress a^) — (a xy )^> /V from Eq. (10). Also, because mostly 
finite wavevectors will be considered, the nonzero averages at q = often will be 
suppressed so that we have for fluctuations l Sfq = / q '. 

Similarly, wavevector-dependent steady-state structure functions from Eq. (4) 
obey: 5/ kffq = NS fg]Cl S k .q, where 

SfUi) = ^ (*/q *5q) W ■ (23) 

The familiar structure factor built with density fluctuations shall be abbreviated by 
"SqCiO = 7r (Qq £q)^- While these findings are familiar from systems without shear, 
translational invariance of sheared systems takes a special form for the two-time 
correlation functions from Eq. (4). Because 

^/ kSq (*)=e- i(q ' Kt+q - k) - a C/ kSq (*) ! 
as follows from Eq. (21), a fluctuation with wavevector q is correlated with a 
fluctuation of k = q(t) with the advected wavevector q(t) = q + q-K t at the later 
time t; only then the exponential in the last equation becomes unity. The advected 
wavevector's y-component increases with time as q y (t) = q y +jt q x , corresponding to 
a decreasing wavelength, which the shear-advected fluctuation exhibits along the y- 
direction. Taking into account this time-dependence of the wavelength of fluctuations, 
a correlation function characterized by a single wavevector can be defined, which 
resembles the equilibrium quantity: C/ kS (i) = A r C/ g;q (t)i5 q ( t ) j k with: 

cw*) - ^ WSw ^ * toq> Ci0 . with = q+ • (24) 

Picking out density fluctuations Qq(t) again, the abbreviation C q (i) = 
jf (^q(t) e ° 1 Qq)^ f° r the intermediate scattering function under shear will be used. 
Similarly for the susceptibilities from Eq. (4) one finds X/k9qW = -^X/s;q(*)^q(t),k 
with the result: 

i 

where the specialisation to density variables shall be denoted by Xq(^) = Xgg--q(^)- 

While these expressions are easily formulated, they suffer from a lack of knowledge 
about ty s . Thus in the following a formal framework is developed within which to 
approximate \& s . 

3. Transient dynamics approach 

The following situation shall be studied: The system is in equilibrium at times t < 0, 
when instantaneously a constant shear rate 7 is turned on: 

°™-{2#:;!!: <-» 
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so that the distribution function at t = coincides with the equilibrium one, 
*(r, t = 0) = * e (r). The solution of Eq. (26) is easily found for t > 0: 

V(T,t) = e n(r )** e (r) . (27) 

The switching on of a real rhcomctcr is supposed to influence the initial variation of 
^>(t) only, which will be neglected in the following as the stationary state, presumably 
reached for t — > oo, will be considered. Rewriting the exponential function, 

e ot = 1 + / dt' e nt ' n , 
Jo 

leads together with Eqs. (5,9) to the formal result for the steady state distribution 
function (where physical units are restored, and the adjoint SO is introduced acting 
on the variables to be averaged with ^ s ): 

* s (r) = tf e (T) + T 1 7 = dt' * e (r) a xy e a 'W . (28) 

This simple result is central to our approach as it connects steady state properties 
to time integrals formed with the shear-dependent dynamics. Knowledge about slow 
relaxation processes in the system can enter. Consequently, the steady-state averages 
from Eq. (22) are given by 

/•OO 

/(7) = (/ q=0 )/y+l J dt' (a xy e nV / q= o) , (29) 
while corresponding expressions hold for the structure functions from Eq. (23), 

• poo 

Sfg-Ai) = (Sf* S 9ci }/N + -L j dt' (a xy e nh ' Sf* Sg^) , (30) 
the fluctuation functions from Eq. (24), 

poo 

<W*) = (Sf* (t) e at 4 SgJ/N + ^ J dt' (a xy e n < «' <5/* (4) e n < « 5 flq ) ,(31) 
and the susceptibilities from Eq. (25) 

X/ fli q(*)= -(Sf* it ) nJe ntt 5 5q )/JV 

• P OO 

- % I dt ' <E ( F * + d *) a *y e *' */q(t) d * ^ ' <W • ( 32 ) 

Note that the averages in Eqs. (29) to (32) can be performed with the known 
equilibrium distribution function. When studying the nonlinear rheology of simple 
fluids, transient correlation functions related to Eq. (29) were found useful in 
thermostatted simulations [15] and in mode coupling approaches [14]. 

3.1. Translational invariance reconsidered 

The time-dependent distribution function \l/(r, t) from Eq. (27) can be used to show 
that a translationally invariant equilibrium distribution function 4' e (r) leads to a 
translationally invariant steady state distribution 4' s (r). To that end, as in Sect. 
2.2.3, ^(r, t) is considered at the shifted positions, r — > T' with r'j = r; + a for all i: 
*(r',t) = e "(r)t-p-K-at ^(p) ^ ^33) 
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where \I/ e (r") = 4 , e (r) was used. The SO fl and the operator P k- a with P from Eq. 
(20) commute, because the shear rate tensor satisfies k ■ k = 0, and because the sum 
of all internal forces vanishes due to Newton's third law: 



(P-K-a) - fi (P-K-a) = 

- [(a-K T ■ K T -dj)] 1=0. (34) 



Therefore, the Baker-Hausdorff theorem simplifies Eq. (33) to 

*(r',t) = e n ^ t e - p - K at * e (r) 

= e n ( r )*e-(S < p 0- K - at * e (r) 

= e n ( r)t tt e (T) , (35) 

where the last equality again holds because the sum of all internal forces vanishes. 
Therefore, 

*(r',t) = *(r,*) (36) 

holds, proving that the time-dependent and consequently the stationary distribution 
function v E , s (r) = lim^oo 4 r (r, t) are translationally invariant even though the SO 
from Eq. (1) itself is not. This applies, at least, in cases without spontaneous 
symmetry breaking. Formally, the role of such symmetry breaking is to discard 
some parts of the steady state distribution function and keep others (with the choice 
dependent on initial conditions). The distributions developed here discard nothing, 
and would therefore average over the disjoint symmetry-related states of a symmetry- 
broken system. 

3.2. Particle conservation 

For the present approach to be tenable, the time-dependent distribution function in 
Eq. (27) needs to approach ^ s at long times. Putting aside at first questions of non- 
ergodic glassy dynamics, the consequences of the conserved variables in the system 
need to be addressed. The particle number is the only conserved quantity, and its 
microscopic conservation law reads from Eq. (19): 

a t e q (t) = nt 0q (t) = i 9 jH( t ) ) (37) 

with the particle flux given by the longitudinal element of the stress tensor (this holds 
generally for overdamped motion as the velocity is proportional to the force) and the 
drift flux, jq(t) = -iq cr q (t) - iv£(t), with: 

crJi(*)=- C nt 'E( 1 -3q-F. 
„«(i) = e nU ]T e* q r < . (38) 

i 

To verify that particle number conservation does not prevent decay of the 
dynamics in Eqs. (28) to (32), the slow density fluctuations are eliminated using 
the equilibrium projection operator [2]: 

P = g< l )(NS q )- 1 (g* (l where PA = g q > (39) 
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with complement Q = 1 — P, where sums over wavevectors q are suppressed because 
of orthogonality. Here, the equilibrium static structure factor is abbreviated as 
Sg = Sqlj = 0) = (f? q £q)/-/V (it will appear repeatedly in the following), and 
idcmpotcncy P 2 = P is clear. 

The correlation functions in Eqs. (29) to (32) can be abbreviated by {<r xy e Q * X), 
with X = f in Eq. (29), X = 8f* <5. 9q in Eq. (30), and X = 5f* (t/) e nV Sg^ in Eq. 
(31); employing the projector P and recalling (o~ xy Qq) = 0, these become: 

(a xy e nH X) = (a xy Q e nH Q X) + (a xy e nH g,) . (40) 

In the first term on the right hand side, already only fluctuations which are 
perpendicular to the hydrodynamic densities appear. The projector Eq. (39) can be 
used to show the vanishing of the second term, where the slow hydrodynamic modes 
could enter. From Eq. (A. 4) in the Appendix (a xy e° f * g q ) = follows, as does 
the fact that reduced dynamics and full dynamics agree for the correlation functions 
needed in Eqs. (29) to (32), viz.: 

(<j xy e at * X) = (o~ xy Q e nh Q X) = (a xy Q e QQ * Qt Q X) . (41) 

The result is perhaps not surprising. The fact that density fluctuations are 
independent of the applied velocity field causes the dynamics leading to the changes 
in the stationary expectation values to be orthogonal to linear density fluctuations. 



3.3. Generalized Green-Kubo relations 

The result Eq. (41) obtained from considering the conserved density fluctuations 
completes our derivation of generalized Green-Kubo relations. While the familiar 
Green-Kubo relations of linear response theory connect transport coefficients to time- 
integrals over projected fluxes [20], Eqs. (29) to (32) enable one to calculate the steady 
state properties of strongly sheared states far from equilibrium. Because of Eq. (41), 
the conserved density fluctuations do not contribute and the dynamics contains no 
hydrodynamically slow mode. 



4. Transient density fluctuations 

The problem of calculating steady state averages is thus converted into one of first 
finding the transient dynamics after switching on the rheometer, and then integrating 
this in order to use Eqs. (29) to (32). The transient density fluctuations will be 
important in this process (see Rcf. [9]) and thus shall be simplified first. Because 
of the equivalence of the particles, the normalized transient collective intermediate 
scattering function can be written as: 

= <4(t) e0t %> = jr K e-*" 6 -"' e nt ^ q ) , (42) 

where p q = e lq l " s is the density of a single tagged particle, which is identical to 
the others. By this trick of singling out a particle, the motion of the surrounding 
particles due to the imprinted flow profile can be specified exactly in the time evolution 
described with one SO. By differentiating one finds: 

d t e -"^' r =' e nh = (-iq-KT a +e -^ K -^ t fit e <q-«T.t) e -<q-K-r.« gfi+t _ 
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Integrating in time, the time evolution operator incorporating advection, in the case 
of density fluctuations, is found as: 

1 r'drU.Cr) 
*q(*) = Q~ (^q* G + 2q> . Where 

fi a (t) = - iqnr s + e - ivK ' T ° t fit e *q-«-r.t t (43) 

and e+ is the time ordered exponential, where earlier times appear on the right. The 
time evolution operator can be worked out explicitly (but to little avail): 

fl s (t) = fi f - iqnr s + iq-K- (2d s +F S ) t- q-K K T qt 2 , (44) 

and because of shear-advection is explicitly time-dependent. 

Projection operator manipulations [14] simplify the time ordered exponentials, 
see the Appendix, and from Eq. (A. 8) follows an exact Zwanzig-Mori type equation 
of motion for the transient density correlators: 

+ r q (t) $ q (t) + / dt' M q (i, $,(0 = . (45) 



Jo 

Here the static projector P s = £j q )(l/S g )(f? q * was employed; because of the equivalence 
of the particles it again satisfies P 2 = P s . The time-dependence of fi s (f ) leads to a 
slightly more general time dependence in Eq. (45) than familiar in equilibrium. The 
'initial decay rate' from Eq. (A. 9) is 

( gq * n s (t) gq ) 

1 qv r J — q 
Oq 

q 2 + q x q y it 2 . 2 2 . as 9 

= +(fc%7tH7*)-^7 ¥ , (46) 

which recovers the 'Taylor dispersion' familiar for non-interacting particles [2]; for 
non-interacting particles M q (f, t') = holds. 

The memory function in Eq. (45) is given by 

M q (t,f) = - (A£(t) U s (t,t') B^(t'))/S q , where 

U 3 (t,t') =e J / 1 ' y , (47) 
and where the generalized longitudinal stress tensor elements are 

(A^(t) = (g q * n s (t)Q s 

B'(t)) =fi s (<)g q ) , (48) 

see Eq. (A. 8). The memory function M q (t, £') encodes the after-effects of the variables 
not treated explicitly in <& q (i) that provide a bath for the density fluctuations. 

In the context of mode-coupling theory, much depends on how this intractable 
object is approximated. There is ample evidence, for dense colloidal dispersions close 
to equilibrium, that a Markovian approximation for M is quite insufficient [21]. Also 
evidence from careful dynamic light scattering tests of mode coupling theory close to 
equilibrium [22, 23] suggest that following Cichocki and Hess [24] a second projection 
step is required. We perform this step now; further discussion is deferred to Sect. ^.1 
below. 

In the second projection operator step, the time evolution operator, Eq. (44), is 
formally [24, 25] decomposed as 

n s (t) = ni(t) + n s (t) gq ) (g q * n s (t) g q )- J <g q * n s (t) 

= ni(t) + n™ d -(t) , (49) 
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where the 'reducible' part of the SO couples the dynamics back to the generalized 
stress elements: 

n™ d -(t) Q s = -B»(t)) — (A s *{t) . (50) 

The importance of this separation of f2 s (t) lies in the possibility to introduce another 
memory function. While Mq(t,t') plays the role of a generalized diffusion kernel, the 
new memory function m q (i, t 1 ) plays the role of a generalized friction kernel. As shown 
in the Appendix, see Eq. (A. 10), the original memory function can be rewritten using 
Eq. (49) as 

M q (M') + r q (t) m q (M') r q (0 

+ r q (t) Jdt" m q (t, t") M q (i", = 0, (51) 

where the new memory function is defined as: 

m q (M') = T~\t) (A°*{t) Ui(t,t') B q (f)> r-^t') /S q . (52) 
Its time-dependence is given by the 'irreducible' [24, 25] dynamics introduced in Eq. 
(49): 



Ui(t,t') = 




From the theory of Volterra integral equations [26], see Eq. (A. 10) in the Appendix, 
it follows that the equation of motion, Eq. (45), thus can be rewritten as: 

d t $ q (t) + r q (t) |# q (t) + J dt' m q (t, t') d t , $ q (i')j = . (54) 

Equation (54), together with the definition of the memory function m in Eq. (52), is 
the central new result of the approach to shear thinning introduced in Ref. [9] , and is 
derived explicitly here for the first time. Together with the generalized Green-Kubo 
relations of Eqs. (29) to (32), it will be the starting point for factorizations building 
on the insights of mode coupling theory [10] into the dynamics of quiescent colloidal 
dispersions. 



4-1- Discussion of the memory functions M and m 

The equations of motion containing the two memory functions differ because of the 
Cichocki-Hess projection step in Eq. (49). The different contents of Eqs. (45,54) can 
be seen from performing a Markovian approximation in the two memory functions. 
Then, Eq. (45) becomes 



r q (*) 



dt' M q (t,i') 



$ q (t) = 



(55) 



where the renormalization of the decay rate T can be expected to be negative (without 
shear this can be shown rigorously, and is connected to Eq. (8)). In order to 
describe a slowing down of the dynamics (viz. a small effective decay rate) thus a 
near cancellation of the two terms in the square bracket in Eq. (55) is required. 
Approximations will need to be subtle to recover this near cancellation. 

On the other hand, Eq. (54) becomes upon performing a Markovian 
approximation: 



fit*q(t) 



r q (t) / i + r q (t) / dt'm^t') 



* q (t) = , (56) 
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where the renormalization of the decay rate now describes a suppression of the 
dynamics (slowing down) as the memory contribution can be expected to be 
positive (again without shear rate this can be shown rigorously), and large. Any 
approximation yielding a large memory-integral (as expected close to equilibrium in 
dense dispersions) , thus can reasonably describe slowing down using Eq. (54) without 
running the risk to predict an unstable system, viz. negative decay rates. 

4-2. Neutral or vorticity direction 

In the vorticity direction, q = q z, perpendicular to the impressed solvent flow and 
its gradient, the Eqs. (52,54) simplify to almost the known ones from the standard 
Zwanzig— Mori approach. The equation of motion becomes 

9t*qt(t) + j- 1 %i{t) + jf dt' m qi (t,t') cV = , (57) 

with the simpler expression of the memory function: 

m qi (t,t') = {(a% +V$") v%)S q /N, where 

v n = J2e^ F J^LQ. (58) 
51 q 

The stress tensor a q was defined in Eq. (38). To simplify m the equivalence of the 
particles was used, replacing the single particle fluctuation by the collective one 
Qq/N in all averages where the index s appears only once. Importantly, shearing 
affects the vorticity direction not only via V K , but also via the reduced dynamics 
which couples all spatial directions and contains the shear rate 7 in any order: 

Besides the reassurance that the formal manipulations recover almost standard results 
in the case where shearing affects the particle motion least, the result Eq. (58) 
for the memory function in the vorticity direction is noteworthy for two reasons: 
First, the stress-stress autocorrelation function m calculated from Eq. (1) without 
shear arises from potential interactions and thus approaches a constant for vanishing 
wavevector, m q (t) — ► mo{t) < 00 for q — > [21]. With shear, however, the result that 
? 2 m q i(t, t') — * const, for q — > can be expected from Eq. (58), because the particles 
are forced by the flow field. Thus the hydrodynamic collective diffusion process will 
be affected. Second, the complicated time dependence which arises in the memory 
function m q (f , t') in x- and y-direction because of the advection of stress fluctuations 
with the imposed flow, simplifies to a dependence solely on the time- difference between 
stress fluctuations along the vorticity direction, m q &(t,t') = m q z(t — t'). 

5. Discussion and outlook 

The derived generalized Green-Kubo relations of Eqs. (29) to (32), and the equations 
of motion for the transient density correlators Eq. (54) follow from the Smoluchowski 
equation for Brownian particles under uniform imposed shear as given in Eq. (1). 
From the Green-Kubo relations, general conclusions about linear response around 
the equilibrium state can be made by setting 7 = in the dynamics. The reduced 
dynamics in Eq. (41) contains no hydrodynamic components, because the density 
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is the only conserved variable in Brownian systems. As long as the equilibrium 
fluctuations are ergodic and decay faster than 1/t for long times, the leading change 
in any stationary variable is linear in shear rate 7. Non-linearities in 7 in steady 
state quantities will be largest for variables where the transient dynamics exhibits the 
slowest algebraic decay 

A central approximation of the approach is hidden in our postulate that the 
time dependent solutions to the Smoluchowski equation approach the stationary 
solution at long times. Aging effects [27] could prevent glassy quiescent states to 
follow the transient dynamics calculated above. Spatial symmetry breaking could 
lead to inhomogeneous states, such as 'shear-banded' ones. Only comparison with 
simulations and experiments can determine whether systems exist exhibiting the 
postulated properties. 

A central difference to standard equilibrium Zwanzig-Mori equations for density 
fluctuations is the appearance of the time-dependent wavevector advection in the time- 
evolution operator fl s (t) of Eq. (43). It arises because of the affine deformation of 
fluctuations and is an exact consequence of shear in the Smoluchowski equation of Eq. 
(1). As found for simple liquids [16, 17], the mode coupling approximations described 
in Refs. [9, 12, 13] deduce from it that shear speeds up the structural relaxation and 
thus causes shear-thinning. The aspect that the stationary nonequilibrium state is 
characterized by a non- vanishing probability current, which is connected to the non- 
Hcrmitian nature of the Smoluchowski operator, enters our approach in the strategy to 
calculate the steady state distribution function via integrating through the transient. 

As discussed in Sect. ^J, the memory function m q (t, 1/) from Eq. (54) and the 
equation of motion Eq. (52) appear reasonable starting points for approximations 
capturing the slow dynamics in driven (sheared) dense colloidal dispersions. Mode 
coupling approximations had been suggested in [9] , their universal contents had been 
discussed in [12], and their detailed presentation will be given in a future companion 
publication. 
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Appendix A. 

The appendix contains various more technical manipulations, which are used in the 
main text. 

Appendix A.l. 

The calculation of the linear response susceptibility in the stationary state starts from 
the change in the energy given in Eq. (7). The SO changes to — Af2(r, t), where 



To linear order in the external field h e the stationary distribution function changes to: 




(A.l) 




(A.2) 
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which leads to the shift of an arbitrary expectation value linear in the external field 
h e given by: 

(g) (A < Jle) - (g) h) = - [ dt' g(T) e ° (*"*') Afi(rt') tt s (r) + <D{h 2 e ) . (A.3) 



A partial integration leads to Eq. (6) with the definition of the susceptibility in Eq. 
(4). 

Appendix A. 2. 

In order to show that conserved density fluctuations do not prevent the dynamics in 
Eqs. (29) to (32) from relaxing, the following operator equality is useful, where the 
first line can be shown straightforwardly by differentiation: 

e uh = e nt Qt + f dt , e nh> nt P e nt Q (t-f) 



= ^ Qt + ^ q ] f eUU ' 3 ^ ^ Q (t "° ' (A ' 4) 

where P is the projection operator from Eq. (39). Two conclusions can be drawn. 
First, because (a xy e n 1 g q ) in Eq. (40) by translational symmetry (see Eq. (21)) can 
be non-vanishing for q = only, and there only the first term in Eq. (A.4) survives, 
it follows that: 

(o- xy e ot * g q ) = (cr xy e nt Qt g q ) 5^ 

= ((Tory Qq) ^q,0 = , 

as can be seen by expanding the exponential and using Q g q = 0. Consequently, Eq. 
(40) simplifies to 

(a xy e nU X) - (a xy Q e nU Q X) . 

Second, reduced dynamics and full dynamics agree for the correlation functions needed 
in Eqs. (29) to (32), viz.: 

(a xy Q e nH QX) = (a xy Q e° + Q 4 Q X) , (A.5) 

again, because of the vanishing of the difference (arising from the second term on the 
last line in Eq. (A.4) ) at q = 0. This result leads to Eq. (41). 

Appendix A.3. 

In order to derive the equation of motion for transient density fluctuations, the time 
ordered product in Eq. (43) is rewritten using the projection operator Q s , which is 
the complement (viz. 1 = P s + Q s ) to P s defined below Eq. (45): 

e{° =U s (t,0)+ dsU s (t, S )n s (s)P s e J + , (A.6) 

Jo 

with the abbreviation U s (t,t') from Eq. (47) . Equality can be shown by differentiation: 
d t e{° = n s (t) Q s U s (t , 0) + n s (t) P s e{° 

+ fds n s (t)Q s U s (t,s) n s (s)P s e{° dTfMT) , (A.7) 
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where d t U s (t,t') = fi s (t) Q s U s (t,t') was used. Regrouping on the right hand side, 
shows that (<9 t — Q s (t)Q s ) rj{t) = Vt s {t) P s cxp + J Q dr fi a (r) , where r/(t) abbreviates 
cither the left or right hand side of Eq. (A. 6). Yet, Eq. (A. 7) turns out more useful 
when sandwiched between density fluctuations: 

r*drO,(r) (nf*n s (t)n a ) pdrrisfT) 

dt $ q (i) S q = d t e{° ; g q ) = ^ q Km/ (4* e{° e,) 

+ /ds (e q *r» s (t)Q s tf,(t, fl ) O s ( s ) eq )i(^ q * e {° dT " s(r) g q ) . (A.8) 

JO Oq 

Here, Q s C/s(£, 0) Qq = Q s Qq = was used. Equation (45) follows and the definitions of 
the rate r q (i) and the memory function M q (i, i'). The rate r q (f) = —S q {g^* fi s (t) £?q) 
can easily be evaluated, using Eq. (44): 

K n,(t) f? q ) = <£ q * (fi e - q-K-Z^-qt 2 ) ft,) 

+ iq-w (g q * (29 s + P.) g q ) + ^<£> q {v V K T -di - q-KT s 5i, s ) g q > 

i 

= - 9 2 -q-K-/ t T -qt 2 5 (; -q^-qi(l + 5 9 )+^ i q^-<f ?q *(r^ s )e 4qr -) , (A.9) 
which leads to the stated result, because the last term in Eq. (A.9) becomes q-K--Jj^S q . 
Appendix A. 4- 

The decomposition of fi s (i) in Eq. (49) leads to the differential equation for the 
reduced dynamics: 

d t u s (t,t') = n s (t)Q s u s (t,t') = n l s (t)Q s u s (t,t') + ni cd (t)Q s u s (t,t') , 

which can be viewed as differential equation with £l r s cd (t)Q s U s (t, t') as inhomogeneity. 
It can be integrated to give 

U s (t,t') = U l s (t,t') 

dt" ui(t,t") B s q(t")) ^T^(t") (A^(t") u s (t",t') , 

where the explicit expression for the reducible part from Eq. (50) was used, and the 
'irreducible' fast dynamics U l s (t,t') corresponds to the solution of the homogeneous 
equation. It is given in Eq. (53). Inserting the expression for U s (t,t') into the 
definition of M q (f , t') in Eq. (47) immediately gives Eq. (51) with the definition of 
the memory function Eq. (52). The equation of motion Eq. (45) can be viewed 
as a Volterra integral equation of second kind for 3> q (i), with kernel proportional to 
Mq(t,t') and — <9 f $ q (i)/r q (t) as inhomogeneity: 

3 q (f) + J o dt' ^Mqff.O S q (f) = d t %(t) . (A.10) 

The solution is given by: 

*q(«) = -fryj - J dt 1 m q (i,i') r q (t') ^-^y a t ,$ q (t') , (A.ll) 
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where the resolvent kernel m q (£, t') satisfies the integral equation [26]: 
M q (i,i') + m q (i, £')r q (i') 



r q (t) 



' 1 



+ Jdt" m q (i,i")r q (i") ^-^y M q (i",i') = , (A.12) 

By comparison of Eqs. (A. 12, 51), the memory function m q (i,i') of Eq. (52) is 
identified to agree with the resolvent kernel entering in Eq. (A. 11), m q (i,i') = 
iriq(t,t'). Thus Eq. (A. 11) actually is equivalent to Eq. (54) as was to be shown. 
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